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in computing energy-momentum tensor correlators on the lattice from which the viscosities can 
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methods should soon yield predictions at the higher temperatures that will be explored at the 
LHC experiments. 
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1. Introduction 



Heavy ion collisions at RHIC have produced a rich phenomenology - see the contribution 
of R. Pisarski at this conference - which has revealed unexpected properties of the quark gluon 
plasma. Here I will be concerned with the thermodynamic and transport properties of the plasma 
in the range of temperatures T c <T < 4T C . Hydrodynamics calculations [jl[ ^ ^] successfully de- 
scribed the distribution of produced particles in heavy ion collisions at RHIC @,|, This early 
agreement between ideal hydrodynamics and experiment has been refined in recent times. On the 
theory side, the dissipative effects of shear viscosity 77 have been included in full 3d hydrodynam- 
ics calculations [||, ||, [T^] and the sensitivity to initial conditions quantitatively estimated JTT| ] for 
the first time. On the experimental side, the elliptic flow observable v%, which is sensitive to the 
value of 77 in units of entropy density s, is now corrected for non-medium-generated two-particle 
correlations [12]. The conclusion that r\/s must be much smaller than unity has so far withstood 
these refinements of heavy-ion phenomenology [pT|]. 

On the theory side, it is therefore important to compute the QCD shear viscosity from first 
principles to complete the picture. Furthermore, since the heavy ion collision program at LHC will 
probe the quark-gluon plasma at temperatures about a factor two higher, it is crucial to predict the 
shear viscosity at « 3T C , and from there to predict the size of elliptic flow, before experimental data 
is available. 

A small shear viscosity is a signature of strong interactions, as I will explain in the next sec- 
tion; strong interactions in turn require non-perturbative computational techniques. In this talk I 
present the lattice calculation of the thermal correlators of the energy-momentum tensor in the Eu- 
clidean SU(3) pure gauge theory, and discuss methods to extract the shear and bulk viscosity from 
them. Computationally, the calculation is challenging enough without the inclusion of dynamical 
quarks, and physically, the thermodynamic properties of the QGP not too close to T c do not depend 
sensitively on the flavor content [13]. In perturbation theory Jl4|], the ratio of shear viscosity to 
entropy density is 0(^? and there is only about 30% difference between the pure gauge theory and 
full QCD [15] at a fixed value of a s (t? /s is smaller in the pure gauge theory). It should be appre- 
ciated that this is not much for a quantity which is infinite in the Stefan-Boltzmann limit, and the 
difference is actually reduced when comparing the pure gauge theory and full QCD at a common 
value of T /T c . 



2. The physical significance of viscosity 

Consider a dilute gas enclosed in a box, Fig. [jj A vertical velocity gradient V/L can be 
established by pulling the upper plate horizontally, at the cost of overcoming a friction force per 
unit surface F/S. The macroscopic interpretation of that force is that neighboring layers of fluid are 
being dragged along. The microscopic mechanism is that molecules passing through a horizontal 
unit surface from above carry more longitudinal momentum than those coming from below. This 
generates a net transfer of momentum in the transverse direction. The longitudinal momentum 
of a molecule that passes through the unit surface is characteristic of the average longitudinal 
momentum in a slice located at a separation of order the mean free path £ m j p from the reference 
slice. The downward flux of molecules is clearly proportional to their density n. However the 
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Figure 1: In the linear regime, the friction force per unit surface felt by a layer of fluid is proportional to the 
transverse velocity gradient. The proportionality coefficient is the shear viscosity: F/S — rjV/L. 



relation between the mean free path and the density, £ m f p = cancels that dependence in the net 
momentum flux (a is the interaction cross-section). The viscosity is then given by T] = where 
(p) is the average momentum. The result is independent of the density n, at leading order in a 
density expansion. 

A finite bulk viscosity has the effect of reducing the pressure of an expanding fluid, as com- 



pared to the adiabatic expansion, P = P equ i\. — £V • v [16]. For a given equation of state, it can 



thus reduce the amount of radial flow in heavy ion collisions. It can also play an important role in 



cosmology [17]. In a dilute gas with two-body interactions, the bulk viscosity is proportional to the 



density of molecules, unlike the shear viscosity, and is therefore much smaller than the latter. 



For further qualitative discussions of viscosity, see [18, |15[ ]. 
3. Methodology 

The set of operators that plays a central role in the calculation of the bulk and shear viscosi- 
ties is the energy-momentum tensor. The Euclidean energy-momentum tensor for SU(/V C ) gauge 
theories reads 

T liV (x) = ll v(x) + ^8 ll v0(x) (3.1) 

fyivC*) = \&iivFpoFp a — F^ a Fy a (3.2) 

d(x) = eg F» a F» a (3.3) 

P(g) = -b g 3 + ..., b = ^. (3.4) 

In contrast with Minkovsky space, To,- = doi is an antihermitian operator. In particular, (r , (x)r , (0)) < 
for x ^ 0, and Pj = if d i xToj(x) is the usual momentum operator, for instance Pj\q) = qj\n) for 
a one-particle state. 

The Euclidean correlators to be computed on the lattice are (Lq = l/T) 

C^ v , pa (x ,q,T) =L 5 Q J d 3 x e^ x (T^ v (0)T pa (xo,x)) . (3.5) 

For diagonal correlators, the spectral functions are then defined by 

s F°° , s coshtofilo -*o) 

C /XViMV (x ,q,r) =L\ \ p A1VjMV (w,q,r) . / m , dco. (3.6) 

Jo sinner 1 
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(no summation over p,v is implied here). The spectral functions are positive, p((0,q)/(0 > 0, 
and odd, p(— ft),q) = — p(ft),q). It is straightfoward to obtain a spectral representation for p, and 
thereby to show [[TJ] that it is related to the Minkovsky-space retarded correlator via p(co,q,T) = 
— ^ImGfi(ft),q,r). The Kubo formulas then relate the low-frequency behavior of G#(ft),q,r) 
with the shear and bulk viscosity, see [^(J for a derivation. The shear and bulk viscosities are given 
by [0] 



71 lim 



Pi 2 ,i 2 (o),o,r) 

ft) 



C(T) = %Y lim Pkk ' ee{( ° AT \ (3.7) 



In Ref. [|22|], I defined the following moments of the spectral function (n = 0, 1, . . .): 

co 2n p{co) d 2n C 



{co 2n ) 



d(D 



sinha>Lo/2 dxl] 1 



xq=Lq/2 



(3.8) 



The latter equality implies that they are directly accessible to lattice calculations. 



4. Discretization of T M v 

In the following sections I present results obtained with the 'clover' discretization of the 



energy-momentum tensor [23, 24]. I will use the anisotropic Wilson gauge actions, see for in- 



stance [25]. Indeed it has long been recognized that such a lattice presents advantages for the 



calculation of thermodynamics [25] and thermal correlation functions, in particular in charmonium 



calculations [|26|]. The downside is that more normalization factors are required. 
The field strength tensor F llv (x) is given by 



Fok 



1 



■(GotW-fiw)) 



%a a a z 

Ki = (Qki(x) - Q, k (x)) . 



(4.1) 
(4.2) 



It is the obvious generalization to the anistropic lattice of the definition found in [|7|]. We define 
the normalized operators as follows, 

a 3 a a t d 00 (x) = -j ReTr |<^ Z dT (g , §o) £(«4a?%(*) 2 ) ~ Zo l Zdo(go, &) £ (a*F u (x) 2 ) } , 



ala z 6{x) = --f^ReTr{^ Z eT (so^o)L(4^^ 

aioga a l k k<l J 



e^ a a T 0Q3 (*) 



J log 
2Z, 



^ReTr 



(a^a T F iF 3 i) + {a^a^F^F^) 



(4.3) 



Here £o is the bare anisotropy, and £ the renormalized one. In total, 6 requires two normalization 
factors, and 8^ v requires five of them. 

To fix the normalization factors of $oo and 6 we proceed as follows. Let L be a physical length. 
Choose an ensemble '1' with Lo = L\ = L 2 = L3 = L and an ensemble '2' with Lq/2 = L\ = L2 = 
L3 = L; then Fq = L 4 ((0qo)2 — ($oo)i) is an RGI quantity. The subtraction is necessary on the 
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anisotropic lattice in order to cancel the mixing with the unit operator. Consider now an ensemble 
'3' where Lq = L2 = L3 = L\/2 = L, and compute F\ = L 4 ( (611)3 — (^n)i)- By symmetry between 
direction 1 and 6, Fq = F\ . This condition fixes the ratio of the two renormalization factors. 

Similarly, requiring the equality of theRGI quantities Go = L 4 ((6) 2 — and Gi = L 4 ((6>) 3 - 
(0)i) fixes the relative size of the two normalization factors entering the trace anomaly operator 6. 

To fix the absolute size of the normalization factors, we equate the clover-discretized confor- 
mality measure (6) with the plaquette discretization, for which the normalization factors are known 



in terms of the lattice beta function [ ]25[ ] and the renormalization of the anisotropy Q28Q. For the 
case of 0oo > we obtain its absolute normalization by calibrating its expectation value to the value of 
the entropy density obtained in [^] from the integral method. 

The normalization factor of is determined by requiring that Co3.o3(*o = 1 /2T,0,T) =s/T . 
This equality can be checked explicitly in hydrodynamics and in free field theory. 

4.1 Cutoff effects in perturbation theory 

The cutoff effects affecting the Euclidean correlators can be studied, and, eventually, removed 
order by order in perturbation theory. 

In the continuum and infinite volume limit, the two-point function of 5(711 — T22) equals the 
two-point function of Tu for all q = (0,0,g). On the infinite lattice, discretizing either form yields 
a correlator that approaches the continuum limit with 0(a 2 ) discretization errors. Figure ||] shows a 
comparison of the discretization errors affecting these two discretization schemes at £ = 1. Using 
the clover discretization of 7^ v at a given value of N t , discretizing T12 yields smaller discretization 
errors than discretizing \{T\\ — T22), for all values of xq. The figure also shows that the cutoff 



effects of the plaquette discretization used in [24] are almost identical to those obtained with the 
clover discretization of T\2- 

In order to find the optimal range of anisotropics, we consider the cutoff effects on the tensor 
correlators at a fixed value of the spatial lattice spacing, <%/Lo =fixed. We then vary the anisotropy 
between 1 and 4. On the right panel of Fig. ^, we see that the sign of the cutoff effects changes 
for 1 < i§ < 2, and goes to a finite value in the Hamiltonian limit, t, —> 00. It is clearly seen that 
any choice E, > 2 reduces the cutoff effects significantly. It also appears that choosing % > 3 does 
not reduce the cutoff effects further, because they are dominated by the coarseness of the spatial 
discretization. The cutoff effects are minimal near % = 2, although, because of the sign change, 
this is partly accidental. 

It should be noted that Eq. ^ holds in the continuum, and therefore the Euclidean correlator 
should in principle first be extrapolated in the standard way, before one attempts to reconstruct the 
spectral function. In practice, we often use data at finite lattice spacing instead, but one must then 
be careful not to use data at too small xq/o, where cutoff effects become large. In the data presented 
in section 6, we have reduced cutoff effects by dividing the non-perturbative lattice correlator by 
the treelevel lattice correlator and multiplying it by the continuum correlator. 

5. Analytic knowledge of the spectral function 

In this section we gather some of the known useful facts about the spectral functions. 
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Figure 2: Left : The relative deviation of the p = tensor correlator at finite lattice spacing from the contin- 
uum correlator. The spatial volume is infinite. The D's are for the clover discretization of T\2 and the (g) for 
the clover discretization of j (T\\ — T?i). The V is the plaquette discretization used in [Bjh. Right: The cutoff 
effects on the tensor correlator for xq —Lq/2. The spatial lattice spacing a a is held fixed, the temporal lattice 
spacing a T is varied between a quarter and one times a a . The D's refer to T\2 and the <g>'s to \(T\\ — T?i)- 
For comparison, the cutoff effect on an isotropic N T = 16 lattice is indicated by a cross. 



5.1 Ward identities 

The Ward identities that result from the conservation of the energy-momentum tensor, d^T^y = 
0, imply, for q = qe^, 

» 4 p 00 ,oo(w,q) = -«Vp03,03(«,q) = <? 4 p33.33(«,q) (5.1) 
-W 2 poi,oiO,q) = q 2 p l313 (co,q). (5.2) 

The first set of spectral functions is referred to as the 'sound channel' and the second as the 'shear 



channel' Q20[]. We therefore use the notation p sn d(w, q) = Poo,oo(&>, q) and p s h(w, q) = Poloi (w, q). 

At zero temperature, Lorentz invariance imposes strong constraints on the analytic form of the 
spectral function. This is shown in some detail in the Appendix for p sn d(w,q). In particular, the 
Euclidean correlator Coo,oo vanishes as |q| 4 at T = 0, whereas it generically only vanishes as q 2 at 
finite temperature. 

5.2 Weak coupling predictions 



In [29], I computed the spectral functions in the free field theory approximation. They corre- 
spond to the result for photons, times the color factor oIa = N 2 — 1. At finite momentum, they are 
functions of two dimensionless parameters, Txq and q/T , and are linear combinations of polylog- 
arithms. An example in the sound channel is shown on Fig. |5[ Unlike in other channels illustrated 
in p9Q, the spectral function has a discontinuity at threshold, a> = q. At zero-momentum, the 
spectral functions simplify considerably. For instance, 



p 33 ,33(«,0,r) = 1g/ A 2 — rw+ TT 3d A T 4 (0S((D) , (5.3) 



2d A ft) 4 (2n\ 2 
15(47r) 2 tanh^ + Vl5/ 
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These free-field theory spectral functions of the fluxes contain an (o8((o) term, which expresses the 
conservation of a momentum flux, an unrealistic feature from the QCD point of view. To describe 
the dissipation of momentum flux, a more careful analysis is required, that was been carried out 
in [ |30| , 16]. These studies show that the delta function becomes a peak of finite width, with a long 



tail in the 7^ channel. 

5.3 Predictions of hydrodynamics 

If the momentum q = (0, 0, q) , the hydrodynamics predictions for the shear and sound channel 

p sh (ft»,q) w, q ^o 17 q 2 



are respectively [20] 



CO 71 CD 2 + (r}q 2 /{e+P)) 2 

Psnd(w,q) m,q^0 T s (E+P)q 4 



(5.5) 
(5.6) 



CO % (CO 2 - v 2 q 2 ) 2 + (T s coq 2 ) 2 ' 

where 

is the sound attenuation length. The sound channel spectral function is illustrated on Fig. 0. 

5.4 Strongly coupled J/ = 4 super- Yang-Mills 

The spectral functions of the jV = 4 SU (N c ) super- Yang-Mills theory were calculated by 



AdS/CFT methods in Q3 1| , [20J at infinite 't Hooft coupling and large number of colors. Those 
spectral functions are very smooth, and exhibit precisely the behavior generically predicted by 
hydrodynamics at low frequencies and momenta, as illustrated by Fig. [| The sound peak in p sn( j 
is the dominant feature of the spectral function at least up to q/lizT = 0.6. This is encouraging for 
lattice calculations, since the data presented in the next section already reach non-zero momenta as 
low as q/2nT = 0.125. 

6. Numerical results 

We simulate the = 2 anisotropic Wilson gauge action with the two-level algorithm as de- 



scribed in [32, 33]. We observe that correlators with a small, non- vanishing momentum are more 
accurate than at q = 0. The typical accuracy achieved on the correlators is 1%, and in some cases 
it is as good as 0.5%. Here we show results for the temperatures 1.27, 1.52, 2.52 and 3.15 in units 
of T c , respectively from lattices 12 • 48 3 , 10 • 20 3 , 10 • 20 3 and 8 • 24 3 . 



6.1 Sound channel 

Due to the Ward identities, the correlators are all determined by the spectral function of the 
energy density. In particular, at fixed xq = I /IT, 

Coo,oo(q,T) = / d(0 . ' , (6.1) 
Jo sinhft)/2r 

-C Q3fi3 {q,T) = -j / dcoco , — — , (6.2) 
q 1 Jo smhft)/2r 
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Figure 3: The SU(N C ) gauge theory sound-channel spectral functions predicted by free theory (top). The 
functional form predicted by hydrodynamics (middle), for v 2 s = A, s = ^s$b and 7T S = i. At the bottom, the 
exact result for 2p / (ticIaT 4 ) in infinite-coupling, large-A^ £ SYM theory, computed by AdS/CFT methods, is 
reproduced from (thex-axis variable is w — 0)/2%T and the curves correspond to q/2nT = 0.3,0.6, 1.0 
and 1.5). 
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Figure 4: The sound-channel Euclidean correlators atxo = 1 /2T. 



Energy -momentum tensor correlators and viscosity 



Harvey Meyer 



1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 



-c -i,oi( x o= 1/ ( 2T )>q) /d A 




T=1 .27 T r 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

( q/ 2ttT f 

Figure 5: A shear-channel Euclidean correlator at xq = 1 /2T. 



. . 1 
C33,33[q,T) = 



o sinhft)/2r 



(6.3) 



Computing Co3,o3 and C33 33 thus provides complementary information to Cbo,oo> without having to 
take finite differences in xq, which necesarily leads to larger statistical and discretization errors. 

Figure ^displays the three correlators of the sound channel, as a function of spatial momentum 
q = qei,. One remarkable feature is that below 2T C , |Co3,03 1 increases with q 2 at small momentum. 
This is unlike the free theory prediction. 

An Ansatz that interpolates between the hydrodynamic behavior at low frequencies and the 
perturbative behavior at high frequencies is 



p S nd(w,<?,r) 



tanh(A) ((0 2 -V 2 q 2 ) 2 + (T s (Oq 2 ) 2 



(l-tanh 2 w/2r) + 



2dAq 4 
15(471) 



tanh 2 w/27\ (6.4) 



Fitting to the correlator Coo.oo^jQjT 1 ) at 1.27T C for q/(^-) = 1, V2, y/3 and 2, I obtain TF S 
0.15(10). This is compatible with my previous estimates [24, ^2|] for rj and £. 



In this (over-)simplified determination, I am relying on the fact that the integral of the sound 
peak multiplied by 1/cosh is of the form c v T 2 + 0(q 2 ), and that the q 2 term has some depen- 
dence on the width of the peak, itself given by the sound attenuation length. However, contributions 
from co = 0(T) are not parametrically suppressed, although at these small momenta they are nu- 
merically very small both at weak coupling and in the strongly coupled SYM theory, see Fig. [3| 

A related point is that at order q 2 , next-to-leading order hydrodynamic corrections to the sound 
peak enter as well, and those are not described in Ansatz (|6T4|). For instance, the sound velocity has 



a momentum dependence at next-to-leading order. At weak coupling [ 16] P33.33 (ft), 0, T)/oo has a 
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Figure 6: Momentum and Euclidean time dependence of the momentum density correlator Coi,oi/^a and 
Q3.03 / d-A in the shear and sound channels, corresponding to open and closed circle respectively. Error bars 
are smaller than the symbols. The temperature is \.21T C . 



peak whose width (Do is related to the inverse of a relaxation time [34]. Since p snt i = (<7/&)) 4 p33,33> 
this effectively makes the sound peak narrower, typically like 



+ 



Further progress will come from simultaneously analyzing the xq and ^-dependence of the 



correlator. 



7. Shear channel 

Similarly to the sound channel, at fixed xo = 1 /2T and with spatial momentum q = qe^, 

-C 02 ,02W,r) = / dco . , (7.1) 

Jo smhft)/2r 

Cn,rz{q,T) = / dcoco . . (7.2) 

q l Jo sinhw/2r 

Figure || displays the momentum correlator of the shear channel, as a function of spatial momen- 
tum. The derivative ^C(^,q = 0,T) is negative, unlike the free theory prediction. 



8. Discussion 



In the two-point functions of the momentum densities Tok, I found that below 2T C the first 
moment, ■£j I Cokfik(^T ,q = 0, T), is positive in the sound-channel (k = 3), and negative in the shear 
channel k = 1. In free field theory, these signs are reversed. There is thus a qualitative difference 
between these non-perturbative correlators and the correlators in the weak coupling limit. By 
contrast, the correlation functions of the fluxes at vanishing spatial momentum are numerically not 
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very different from the free-field prediction [^4|, 22, 19^. This is because they are dominated by the 
perturbative ft) 4 term in the spectral function at high frequency. 

I pointed out the advantage of the conserved-charge correlators in terms of sensitivity to the 
low-frequency region co <T region [29]. More generally, it is advantageous to constrain simulta- 
neously the ft) and q dependence of the spectral functions, since hydrodynamics makes a combined 
prediction in these two variables. Figure |6] displays the simultaneous xq and q dependence of the 
momentum density correlators in the shear and sound channels. Thanks in part to treelevel im- 
provement, the q = correlator is x$ -independent to very good accuracy, as it should be. It is by 
fitting the whole surface parametrized by (xo,q) that one can best exploit the lattice data. 

Recently, I have studied spatial correlation functions of the energy-momentum tensor [|35|]. 
The advantage of these observables is that no analytic continuation is required: they directly contain 
information on the space-like correlations present in the quark-gluon plasma. They can also be 
related to the spectral functions: 

C(r,r) = lim r d coe-™ P ^l (8.1) 

v ' e^oJ (2n) 3 Jo tanhft)/2r' 

which illustrates the wealth of information encoded in p(ft),q,r). A hallmark of ordinary liquids 
is the short-range order observed (typically by neutron scattering) in the particle pair correlation 
function, which displays several gradually damped oscillations. One of the conclusions I reached 
is that the fluids with the smallest shear viscosity are those that are dense, but that are not ordered 
too strongly. The thermal contribution to the energy density correlations exhibits a non-monotonic 
behavior on distance scales of order l/T that is markedly stronger than in leading order perturbation 
theory. This may be the signature of a spatial ordering that is the quantum relativistic analogue of 
the short-range order characteristic of classical non-relativistic fluids. An attractive strategy to 
test the claimed liquid-like properties of the quark-gluon plasma [ p^ ] is to compare these spatial 
correlators with those computed in the strongly coupled SYM theory, which is known (from its 
spectral functions and small viscosity) to be an excellent fluid. 

In conclusion, there has been much technical progress since the first calculation of the Eu- 
clidean correlators relevant to shear and bulk viscosity [21, 37, 24, The first accurate data in 
the shear channel allowed me to derive an upper bound on the shear viscosity, T] /s < 1.0 at T < 2T C 
under only weak assumptions [24]. To improve on this, the current emphasis is on exploiting si- 
multaneously all the relevant lattice data, as well as the analytic properties of the spectral function 
known from perturbation theory in the ultraviolet and from hydrodynamics in the infrared. Spa- 
tial correlators provide complementary information that may help us understand the ability of the 
strongly coupled QGP to flow. 
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A. The vacuum energy density correlator 

At zero temperature, we can exploit Lorentz invariance to constrain the analytic form of the 



spectral function. Following [38], 



/ 



d A xe iqx (T^(x)T pa (0)) = n tiV ,pa{q)T{q-q)+f^, a p{q)T s {q-q) (A.l) 

Tljiv = q^qv ~ q 5^v> fnv,po = TlnvTlpo (A. 2) 

2 

^jxv,pa — T]/ipT]v,p + T]/jfjT]vp — ^ TlnvTlpo- (A. 3) 

This implies the following form of the energy density correlator in the mixed representation: 

Coo,oo(*o,q) = H A fJ^e- i ^T{ql + q 2 ) + T s {ql + ^)]. (A.4) 
In particular, the correlator satisfies 

f 37T - A) c oo,oo(*o,q) = | q |*o Coo.oo 0o,q)- (A.5) 
dx Q \d\q\ \q\J 

On the other hand, if 

v n (q) = (n,q\L- 3 / 2 f ^x^ x f 00 (x)|£l), (A.6) 
the Kallen-Lehmann representation of the correlator reads 

Coo,oo(*o,q) = J £\v n \ 2 (q)e- E "^ X0 . (A.l) 

n 

On then easily shows that V«, 

E 2 (q)=E 2 (0)+q 2 , |v n | 2 (q) = ^ ) f' q|4 ; (A.8) 

with a n a dimensionless number. In particular, all the T = contributions vanish as ~ |q| 4 when 
q^0. By the Ward identities of the energy-momentum tensor, C^^xq, \q\) = |q| -4 4*0)0 oo(*o>q) 



Noticing that for a scalar state S, {S^l 2 J d 3 xf 33 (x)\a) = \(S\L- 3 / 2 J d 3 xd(x)\Q.}, we then find 
the matrix element for the lightest glueball state, 

s 

a = T , (A.9) 

3M 3 ' 



where s, the matrix element of the trace anomaly, is defined and calculated in [39]. Similarly, 
realizing that C 3 j :3 j(xo,0,T) = ^Ckk,u(xo,Q,T) + |Ci2,i2(*o,0,r), we find that the matrix element 



d2 corresponding to the tensor glueball is related to the matrix element t defined in [ ]39| ] by 

« 2 = 2V3^. (A.10) 
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